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Abstract 

Spatial clustering and intermittency in the relative velocity of heavy particles of the same size 
settling in turbulent flows can be strongly affected by gravity. We present a model for the timescale 
of the fluid velocity gradient seen by particle pairs and propose an effective Kubo number based on 
this timescale to explain the mechanism of gravity-enhanced clustering. We explore the mechanisms 
of the gravity-induced reduction or enhancement of the intermittency in the particle radial relative 
velocity (RRV) at different Stokes numbers based on backward-in-time relative dispersion and 
preferential sampling of the fluid field. These effects of gravity on clustering and the RRV must be 
parameterized in the geometric collision kernel. 
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Collisions of heavy particles induced by turbulence are a general and crucial phenomenon 


with many applications 


3 , 


and fluidization 


a 


Elfli. 


including planetesimal formation rain droplet formation 


of particle collision 


gravity 


3 QJiel. 
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Turbulent clustering and relative velocity are two important aspects 


12| . Particles suspended in a turbulent flow typically move under 


However, the effects of grav ity on clustering and the relative motion 


have seldom been studied until very recently 


17H21|: thus, an investigation of the effects of 


settling effects on clustering and the intermittency in the redial relative velocity (RRV) [22|] 
is the primary topic of this letter. 

The clustering of heavy particles in turbulent or random flows is caused by centrifugal 


force 


23N25| or multiplicative amplification 2a, l27| , respectively. Particles at a small Stokes 
nged out of regions of high vorticity toward regions of high strain rate 


number are centri 

and low vorticity 23|. By contrast, clustering in a random flow at a small Kubo number 
is determined by the multiplication of many independent random exp ansion or contraction 


factors of small volumes spanned by a triad of particle separations 


27| . If this multiplication 


decreases after a long time, then the fluctuations in particle number density held will be 
amplihed, and thus clusters can be observed. Here, the Stokes number is Stx = Tp/rx, 
where Tp is the particle relaxation time and tk is the Kolmogorov timescale. The Kubo 


number is Ku = u^Tf/f, where f is the smallest characteristic length scale, Tf is t 


le smallest 


characteristic timescale, and uo is the characteristic velocity scale of the how 28 


,l29L 


Previous studies indicate that gravity reduces clustering of particles at small Stx 
whereas, very recent results show that gravity may enhance clustering at large Stx 


30 
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33| . Gravity-driven enhancement of clustering is attributed to a Gaussian delta-correlated 


(in time) huid velocity gradient held or a how held with a small Kubo number; and the rate 


of caustic formation is found to be reduced by gravity 


18 


19| . Here, our aims are to study 


Tg, the timescale of the Lagrangian correlation of huid velocity gradient seen by particle 
pairs under gravity that ahects clustering, and furthermore, to explore the mechanism of 
the remarkable changes in the intermittency in RRV due to gravity. 


The direct numerical simulation plus poin 


particle model is typically used to solve for the 


motion of heavy particles under gravity 34j-l36| . The statistical quantities for the hows are 
listed in Table HI A total of 1.2 x 10® particles are tracked. The dynamics of the particles are 
controlled by the Taylor microscale Reynolds number, Rex; Stx; and the Fronde number, 
Frg = ax/g, where ax is the Kolmogorov acceleration and g is the acceleration due to 
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TABLE I. Statistical parameters in isotropic turbulent flows. is the grid number in the flow 
domain of (27r)^. u' is the root mean square (rms) of the velocity, e is the energy dissipation 
rate per unit mass. Lj is the longitudinal integral length scale, vk, V and tk = are the 

Kolmogorov velocity scale, length scale and time scales, respectively. Kmax = A'/S is the cutoff 
wave number in the pseudospectral method. 


Case 

Ar3 

Rex 

e 

ly 

u' 

i/ 

Vk 

tk 

V^max 

I 

128^ 

73.9 

3434.7 

0.095 

19.10 

0.91 

4.31 

0.0051 

1.10 

II 

256^ 

100.9 

3468.0 

0.049 

19.52 

0.99 

3.62 

0.0037 

1.15 
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FIG. 1. (Color online) The degree of clustering as functions of the Fronde number and the Stokes 
number: (a) The RDF, g{R), at a separation R = O.bry versus the inverse of the Fronde number at 
different Stokes numbers, (b) The correlation dimension Do versus the Stokes number at different 
Fronde numbers. The results reported by Bee et al.(2014) 1^ for Rex = 460, represented by lines 
with open symbols, are plotted for comparison. 

gravity. Because the results for Case I and Case II are similar, we will present the results 
only for Case II. 

Movies showing the effects of gravity on clustering are provided in the Supplemental 


Material 


37j . in which gravity suppresses clustering at Stx = 0.6, whereas it enhances 


clustering at Stx = 5.0. Similar observations can be 
20|. We use the radial distribution function (RDF), g{r) 


bund in recent publications 18 


Do 


11| and the correlation dimension 
18| as measures of clustering. Figure [T] illustrates g{R) {R = O.Br/) and D 2 as functions 
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FIG. 2. (Color online) The timescale of the Lagrangian correlation of the velocity gradient Tg as a 
function of Fr~^ at Stx = 5. Tg decreases and it indeed becomes less than tk as Fr~^ increases. 
The solid line represents a model of the form Tg = [(I/t^^o)^ + 


of the Fronde number and Stokes number, respectively. Our results about D 2 versus SIk 


display similar behaviors with those from Ref. jl8|, as shown in Fig. 1(b) In Fig. 1(a) 


for SIk < 1.0 (Lines (1) — (3)), the RDFs monotonically decrease with an increase in 
Fr~^. For Stx > 1-0, (Lines (4) — (8)), the RDFs initially increase to a maximum and 
then decrease with increasing Fr~^. Therefore, gravity enhances clustering at Stx > 1 and 


reduces clustering at Stx < 1- The vertical dashed line in Fig. 1(a) indicates that gravity 
becomes relevant at Fr~^ > 1. 

The gravity-driven reduction of the clustering at small Stx is due to the ineffectiveness of 
the centrifugal mechanism because the particles do not have sufficient time to be centrifuged 
out of regions of high vorticity in the flow when the particles are rapidly going through an 
eddy under gravity. At large Stx, the reduction in Tg due to gravity causes the gravity-driven 
enhancement of clustering, as demonstrated as follows. The linearized equations for the 
relative motion of a nearby particle pair, denoted by the relative position 6xpi and the relative 
velocity Svpi, are Sxpi = 6vpi, and 6vpi = [Fij(t)6xpj — Svpi\/Tp. Here, the dots denote the time 
derivative d/dt] Fij{t) = dui/dxj{x.^\t),t) is the gradient of the fluid velocity experienced 
by the reference Particle 1; and i,j = 1,2,3, denoting the 3 direction components. In 
Fig. [2], we show how Tg varies with Fr~^ at St^ = 5.0, where Tg = Rgradif)dt = 
ir /{Fij{0)F„iri{0))dt with i = j = m = n = 2 along the direction of gravity. 
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FIG. 3. (Color online) The standardized PDF of the RRVs of heavy particles with a separation 
R = 0.5r/ versus the Froude number:(a) Stx = 0.6 and (b) SIk = 5.0. The insets illustrate the 
skewness of the PDF as a function of the Froude number. 


Note that F 22 (t) = du^./dxo = —{dui/dx^ + du^/dx^), denoting the horizontal convergence 


in an incompressible flow 


20|. The Lagrangian correlation Rgradit) versus Fvg is plotted 


in the inset. This plot shows that Rgradit) and, thus the timescale Tg decay rapidly with 
decreasing Fvg. We observe that Tg can even be less than tk when Fvg is sufficiently small. 
Tg can be approximated using a model of the form Tg = [(1/r,,o)^ + represented 

by the solid line, where I is the length scale of the fluid velocity gradient, I = 5.6r] Oiv), 
and Tg^o is the timescale of Tg when Fvg = cx3. Tg^ = 2Atk ~ OiTx)- Thus, we approximate 
Tg using [(I/tx)^ + = (h/'n’o)[l + ~ h/'^’o when vk -C Wq. Using 

rj, Vk and Tg as the characteristic lengthscale, velocity scale and timescale, respectively, we 
dehne an effective Kubo number for a turbulent flow with heavy particles settling within 
it as Kue = VKTg/rj = Frg/StK- Therefore, KUf. -C 1 for heavy particles at small Fvg 
and large Six and turbulent flow held seen by settling particles behaves as a random how 
with a very short timescale. Multiplicative amplihcation begins to play a role in clustering 
enhancement through many independent random accelerations |^. 

We now turn to the ehects of gravity on the statistics of the RRV. For a particle pair 
consisting of Particle 1 and Particle 2 with a separation R = |R| = |xp^^ — Xp^^|, their RRV, 


WriR) = (V; 


- 


d2)^ 


R, is intrinsically related to the clustering. Here, the superscript 


(1) and (2) denote particles in a pair. R = R/|R| is the unit vector pointing from Par¬ 
ticle 1 to Particle 2. When the level of clustering approaches a statistically steady state, 
the negative skewness, S = denotes the clustering tendency, where 
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FIG. 4. (Color online) BIT relative dispersion of heavy particles as a function of Fronde number; 
(a) StK = 0.6 and (b) Stx = 5.0. 


w'^{R) = Wr{R)/\/{wr{Ry). At Stx = 0.6, as shown in Fig. 3(a), the PDFs of the RRV 
are negatively skewed and the absolnte valne 151 decreases with decreasing Fvg as shown 
in the inset, which means that the degree of clustering also decreases with decreasing Fvg. 


By contrast, at Stx = 5.0, as shown in Fig. 3(b), |5| initially increases with decreasing 
Fvg down to Fvg = 0.197, implying an enhanced degree of clustering. When Frg decreases 
further, |5| begins to decrease, implying that the degree of clustering begins to weaken. The 
variations in the skewness at small and large Stx are consistent with the variations in g{R) 


versus Fvg, as shown in Fig. 1(a) 


Another interesting observation in Fig. [3] is that the tails of the standardized PDFs of the 
RRV exhibit remarkable and opposite changes with decreasing Fvg at small and large StK- 
The tails become increasingly narrower with decreasing Fvg at Stx = 0.6, whereas they 
become increasingly broader with decreasing Fvg at Stx = 5.0. These observations reveal 
that the intermittency in the RRV is reduced at Stx = 0.6, but enhanced at Stx = 5.0 
with a decrease in Fvg. Moreover, if gravity is absent, the PDF tail exhibits a shape of 
the form exp[—G at a large Stokes number of Stx = 5.0, which is consistent with 

the prediction of the variable-range projection model proposed for the relative velocities of 


heavy particles at large Stokes numbers (the long-dashed curve in Fig. 3(b)) |38|. We shall 
interpret this interesting phenomenon from the perspective of the backward-in-time (BIT) 
relative dispersion. 

The RRV of a particle pair with a separation R can be expressed as Wr{R) = t“^R ■ 
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Au(xp^^(t), Il{t),t)exp{t/Tp)dt, where Au(xp^^(t), R(t), t) is the fluid velocity difference 
for a separation vector R(t) experienced by the particle pair at t, Xp^^(t) is the position 
of the reference particle. Fig. 0] shows the rms of the separation increment = 

(|R(t) — R(0)p)°'^ normalized with respect to r], as a function of Frg at Stx = 0.6 and 
StK = 5.0. At Stx = 0.6 and in the absence of gravity, increases within the 

interval [—3rp,0]; however, it remains less than 3t] and thus it is located in the viscous 
subrange. Gravity reduces due to the reduced timescale between particles and 

eddies. The intermittency of the flow seen by the particle pairs changes very little with 
a decreasing because it is approximately saturated at such small separations for 

a given Reynolds number [39[. At Fvg = oo, particles at such small Stokes numbers are 
preferentially concentrated in regions of high strain rate with higher intermittency. Gravity 
pulls these particles out of such highly intermittent regions, and the reference particles Xp^^ 
become distributed increasingly uniformly throughout the entire held, and the intermittency 
of the velocity difference seen by the particle pairs decreases, thus, the intermittency of the 


particle RRV also decreases, as shown in Fig. 3(a) By contrast, at Stx = 5.0 and in the 
absence of gravity, can break free of the viscous subrange and return to the inertial 

and even the energy-containing subranges within the interval [—3rp, 0]. It can be as large as 


80?7, such that the Gaussian assumption regarding t 
used in the variable-range projection model is valid 


le particle relative velocity at earlier time 
. For particles at large Stokes number. 


gravitational settling significantly reduces bringing it from the inertial subrange 

into the viscous subrange. At Stx = 5.0, the clustering and Xp^^ are ergodic 27|. The 
intermittency of the flow seen by particle pairs with a small separation increases. 


as does the resulting intermittency of the particle RRV, as shown in Fig. 3(b) The balance 
between the mitigation of preferential sampling small-scale flow structure at small St^ and 
the great reduction in the BIT relative separations at large St^ results in a transition of the 
variations in the tails of the PDFs occurring at SIk ~ 3.0, according to our DNS data. 

The reduction in BIT relative separation also results in a reduced tCr.rms- Fig. |5] shows 




r,rms 


as a function of Stx at different Fvg. In the absence of gravity, Wr,rms increases 


with increasing Stx, and then slowly decays after reaching a maximum value. At very 
large Stx in the absence of gravity, Wr^rms oc 0, where T^p is the 

timescale of the Lagrangian correlation of fluid velocity experienced by a single-particle. For 
a given Fvg, the curve turns downward following a scaling law of the form rather than 
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FIG. 5. (Color online) Wr^rms as a function of Stokes number at different Proude numbers. At very 
small Stokes numbers, the numerical results are consistent with the prediction of the Saffman and 
Turner limit ^|. 


when settling velocity wq = gTp becomes large enough. This can be interpreted as 
follows: for particles at rapid settling velocities, the timescale is Tip ~ Lf/gTp. Therefore, 

'^r,rms ~ y/Tj /g/Tp ~ ' 

In summary, clustering is monotonically reduced at St^ < 1 under gravity because of 
the ineffectiveness of the centrifugal mechanism, whereas clustering is non-monotonically 
enhanced at Stx > 1 because of the multiplicative amplihcation mechanism in the case that 
the proposed effective Kubo number is much less than 1. The mechanism under the opposite 
variations in the tails of the PDFs of the RRV is explored from the perspective of BIT 
relative separation and preferential sampling of small-scale flow structures. The rms of RRV 
is greatly suppressed due to gravity and it follows a scaling law of the form St~^ at a given 
Frg. This understanding of the remarkable changes in the PDF of the RRV challenges the 
current intuition that gravity predominantly affects the relative velocity between particles 
of different sizes falling at different velocities. Therefore, the effects of gravity on clustering 
and the RRV are a signihcant influence that must be considered in the parameterization of 
the geometric collision kernels of monodisperse particles. 
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